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Abstract 

This  document  describes  the  development  of  the  Hydrodynamic  Power  Semiconductor  Device 
Simulator  (HPSDS),  which  solves  the  quantum  mechanical  hydrodynamic  balance  equations  for 
the  electrons  and  holes,  along  with  the  Poisson  equation,  for  a  two-dimensional  device.  The 
program  will  compute  the  electrostatic  potential,  electron  and  hole  densities,  recombination  rate, 
electron  and  hole  velocities,  current  densities,  electron  and  hole  temperature,  and  other  quantities 
of  interest  as  functions  of  applied  bias.  The  program  allows  extensive  analysis  of  unipolar, 
bipolar,  and  heterostructure  devices.  Power  devices  such  as  thyristor  and  other  four  or  five 
terminal  devices  are  treated  under  all  possible  operating  conditions.  Both  steady  state  and 
transient  operating  conditions  can  be  considered.  Because  of  the  inclusion  of  energy  balance 
equations  and  heat  diffusion  equation,  accurate  high-field  characteristics  can  be  obtained  with 
this  program. 

Silicon  and  Gallium  Arsenide  are  the  primary  materials  considered  in  this  program,  although 
other  materials  can  be  included  if  the  user  supply  the  necessary  material  data  file.  In  particular, 
wide  band-gap  materials  such  as  Silicon  Carbide  and  Gallium  Nitride  are  being  currently  added 
to  the  materials  library  of  the  code.  All  materials  related  parameters  are  contained  within  a  single 
subroutine,  so  that  it  is  easily  modified  by  the  user  to  suit  his  own  situation.  For  the  theoretical 
basis  and  numerical  implementation  of  the  code,  the  reader  is  directed  to  the  references. 
Materials  parameters  are  also  discussed  in  the  references. 


1.  Introduction 


1.1  Scope 

This  document  reports  on  the  project  entitled  “Power  Semiconductor  Simulation”.  The  following 
paragraphs  will  identify  a  software  system  for  accurate  and  computationally  efficient  codes  for 
simulating  high-power  semiconductor  devices.  The  codes  were  developed  for  use  in  any  parallel- 
processing  environment,  which  has  the  Massage  Passing  Interface  (MPI),  a  de-facto  standard  of 
parallel  programming  model. 


1.2  Objectives 


The  goal  of  this  project  was  to  develop  a  working  simulation  package  for  power  semiconductor 
devices,  such  as  large  area  thyristors,  and  high  frequency  devices,  such  as  MESFETs,  HEMTs 
for  use  in  oscillators.  Aside  from  the  traditional  materials  to  be  simulated  such  as  silicon  and 
gallium  arsenide,  wide  band  gap  semiconductors  like  silicon  carbide  will  be  simulated  because  of 
its  good  thermal  and  high  voltage  properties.  The  simulation  code  will  be  based  on  scaleable 
hardware  and  will  be  rapidly  accessible  to  DoD  researchers. 


The  project  of  large-scale  numerical  simulation  of  semiconductor  devices  involves 
computational  tasks  that  are  too  demanding  for  single  stand-alone  computers  and  workstations. 
These  tasks  are  best  tackled  with  massively  parallel  computers.  The  numerical  codes  developed 
in  this  project  will  be  specially  tailored  to  suit  massively  parallel  machines. 


The  code  will  enable  users  to  perform  full  electrical  and  thermal  analysis  of  semiconductor 
devices  under  high  power/voltage/current  operating  conditions,  thus  making  it  a  very  useful 
design  and  simulation  tool  for  power  semiconductor  applications  as  well  as  applicable  to  self¬ 
heating  device  studies  and  temperature  dependent  model  generation. 


1.3  Approach 

Modeling  of  power  semiconductor  devices  is  more  challenging  than  the  modeling  of  ordinary 
semiconductor  devices  in  that  the  former  involves  not  only  charge  transport  but  also  thermal 
transport  in  the  devices.  Further,  the  conduction  of  electricity  is  coupled  with  the  conduction  of 
heat.  Thus  the  modeling/simulation  of  power  devices  depends  upon  the  solution  of  a  set  of 
coupled  nonlinear  partial  differential  equations,  the  Poisson  equation,  the  current  continuity 
equations,  and  the  semiconductor  transport  equations,  as  well  as  the  heat  diffusion  equation, 
subject  to  the  appropriate  initial  and  boundary, conditions.  The  task  of  solving  such  a  complicated 
system  of  equations  entails  extremely  high  demand  for  computing  resources,  both  in  terms  of 
CPU  cycles  and  memory.  Massively  parallel  hardware  and  corresponding  software  will  be 
needed  to  accomplish  such  a  task. 


1.4  Major  Program  Elements 


The  package  contains  several  functional  blocks,  each  of  which  is  designed  to  perform  a  certain 
task.  Each  block  is  made  of  many  program  units,  which  perform  sub-tasks.  The  program  units 
have  one  or  more  subroutines  or  functions.  In  addition,  a  library  consisting  of  commonly  used 
tools  are  attached,  along  with  several  modeling  drivers.  The  code  is  written  in  FORTRAN-90 
language. 

The  blocks  are: 

Control  (CF), 

Definitions  (DF), 

Discretization  (DS), 

Devices  (DV), 

Local  Quantities  (LC), 

Materials  Parameters  (MT). 

The  following  is  a  complete  list  of  the  program  units  in  each  functional  block.  All  the  units  are 
coded  in  FORTRAN-90  at  the  present  time. 

The  package  consists  of  several  functional  blocks,  each  of  which  is  designed  to  perform  a  certain 
task.  Each  block  is  made  of  many  program  units,  which  perform  sub-task.  The  program  units 
have  one  or  more  subroutines  or  functions.  In  addition,  a  library  of  commonly  used  tools  is 
attached,  along  with  several  modeling  drivers. 

The  blocks  are: 

Control  (CT); 

Definitions  (DF); 

Discretization  (DS); 

Devices  (DV); 

Local  Quantities  (LC); 

Materials  Parameters  (MT). 

The  following  is  a  complete  list  of  the  program  units  in  each  functional  block.  All  the  units  are 
coded  in  FORTRAN  90  at  the  present  time. 

CT_ADT 

CTDF 

C^JENERG 

CT_PARTI 

CTPOISS 

DFENERG 

DF_MOMEN 

DFPARTI 

DFJPOISS 

DS  CARRIER  DENSITY 


DSCARRIERTEMPR 

DSCHEM3CALP0TENTIAL 

DSCURRENTDENSITY 

DSDOPENDENSITY 

DS_ENERGY_DENSITY 

DSJENERGYFLOWJDENSITY 

DSENERGYRATE 

DS_EQUI_EQS 

DSINVERSEMASS 

DSNPARA 

DSNPARACURRENT 

DSNPARAPOTENTIAL 

DS_NPARA_TEMPR 

DSPARTICLERATE 

DSPHASEVELOCITY 

DS_POTENTIAL 

DS_SCAT_RATE 

DSVELOCITYRATE 

DS_VELOCITY_TENSOR 

DV1D 

DV_2D 

DVBOUND 

DVMyBOLISTICID 

GRID 

HY  SOLVERUNITS 
HyBE 

LCEQJUNITS 

LC_C  ARRIERDEN  SHY 

LCCHEMPOT 

LC_CORR_FUNC 

LC_CURRENT_DENSITY 

LC_ENERGY_DENSITY 

LC_ENERGY_FLOW_DENSITY 

LCENERGYRATE 

LC_EQUI_EQS 

LC_INVERSE_MASS 

LCNPARA 

LC_NPARA_CURRENT 

LCNPARAPOTENTIAL 

LCNPARATEMPR 

LCPARTICLERATE 

LCPHASEVELOCITY 

LC_SCAT_RATE 

LC_SINGLE_LT 

LCSINGLELTIMP 

LC  SINGLE  LT  PH 


LCVELOCITYRATE 

LCVELOCITYTENSOR 

MT_ADT 

MT_A1GAAS 

MT_GENER 

MTJCANE 

MTPARAB 

MT_USER 

WKSPACE 

LIBELEMFUN 

UBJEJPOLATTON 

LIBINTG 

LB_INTG_GAUSS 

LEBJMATRIX 

LB_N_CONST 

LIB_P_CONST 

LIB_ROOT 

LIB  VECP 


2.  Major  Tasks 

The  following  is  a  list  of  major  tasks  with  the  above  observation  in  mind. 

1.  Parallelization  of  existing  code  for  ballistic  diode. 

2.  Single  particle,  time-independent  code  implemented 

3.  Modification  of  code  to  p-n  diode  and  then  power  diode. 

4.  Single  particle,  time  independent  code  parallelized 

5.  Single  particle,  time  dependent  code  implemented 

6.  Modification  of  code  for  multilayer  devices 

7.  Many-body,  time  independent  code  implemented 

8.  Single  particle,  time  dependent  code  parallelized 

9.  Interface  code  development 

10.  Many-body,  time  independent  code  parallelized 

11.  Many-body,  time-dependent  code  implemented 

12.  Single  particle,  time  dependent,  full  band  structure  code  implemented 

13.  Many-body,  time  dependent  code  parallelized. 

3.  Major  Accomplishments 


Schedule  of  Events 

The  following  schedule  reflects  the  schedule  and  status  of  the  major  tasks  outlined  above. 


Module/Tasks 

Expected  Completion 

Actual  Completion 

task  1 

1/97 

1/97 

task  2 

3/97 

3/97 

task  3 

7/97 

7/97  ' 

5.  Alpha  test  summary 

The  following  is  a  summary  of  the  milestones  and  metrics  demonstrated  at  the  off-site  CEN-2 
Alpha  testing  conducted  at  Wright-Patterson  AFB,  on  April  14,  1998. 

•  63%  speed-up  achieved  with  16  processors. 

•  MPI  used  on  SGI  Origin  2000  and  SGI-PCA  supercomputers. 

•  Now  in  process  of  moving  parallel  code  to  IBM  SP-2. 

•  Reasonable  load  balance  achieved. 


Objective 

Required 

Observed 

Scaled  Speedup  (16  nodes) 

40% 

63% 

%  code  in  high  level 
language 

80% 

100  % 

Number  of  HPC  platforms 
software  will  run  on. 

2 

2 

#  of  validations  per  week 

4  or  more 

4 

User  computer  expertise 
required.(l-10) 

Required  =  2  (no  GUI) 
1=GUI,  5=  need  to  modify 
source  code  each  run. 

2 

Major  error  &  %  fixable  per 
week. 

10,  80  % 

5, 100% 

%  Of  software  allows  for 
interactive  processing. 

User  specifies  input  and 
waits  for  output  in  a  non 
real-time  mode. 

User  specifies  input  and 
waits  for  output  in  a  non 
real-time  mode. 

#  of  validations  on  full  code 

1  per  2  months 

1  per  2  months 

%  of  required  external  code 
interfaces  completed 

Greater  or  equal  to  60% 

60% 

%  Of  last  validations  that 
were  accurate  and  valid. 

60% 

100  % 

#  Of  potential  users  exposed 
to  new  software  capability 
through  problems  solved  for 

Greater  or  equal  to  10 

10 

them. 


The  parameters  that  can  be  calculated  by  computation  method  are  summarized  in  the  next  table: 


Carrier  density 

Yes 

Carrier  temperature 

Yes 

Electrostatic  potential 

Yes 

Electric  field 

Yes 

Current  and  current  density 

Yes 

All  outputs  currently  viewed  with  third  party  graphics  visualization  software. 


6.  Beta  Test  Results 


Objective  (CTP) 


Scaled  Speedup 


Scaled  Speedup 


%  code  uses  standard 
higher  order  languages,  etc 


Number  of  HPC  platforms  software  will  run  on  with 
valid  results. 


4  or  more  validation  testes  on  full  code/month 


Levels  of  users  expertise  software  provides. 


#  major  errors  or  <  each  week  and  %  fixable 


#  data/sets/user  scenarios  tested/week 


Beta  Test 


Reduces  clock  time  16x  or  >  on  32 


2, 100%  within  2  workdays 


%  of  possible  useful  output  parameters  provided 


90% 


%  of  required  external  code  interfaces  completed 

REQ 

80%  or  greater 

OBS 

%  of  last  10  software  runs  produced  accurate,  valid 
output 

REQ 

80% 

OBS 

#  of  new  users  added  to  CTA’s  HPC  community  over 
last  year. 

REQ 

15  or  more 

OBS 

References 

1  P.  Zhao,  H.L.  Cui,  D.  L.  Woolard,  K.L.  Jensen  and  F.  A.  Buot,  “Equivalent  Circuit  Parameters  of  Resonant 
Tunneling  Diodes  Extracted  From  Self-Consistent  Wigner-Poisson  Simulation”,  IEEE  Transaction  of 
Electron  Devices,  48, 614-627  (2001). 

2.  p  Zhao,  H.L.  Cui,  D.L.  Woolard,  and  FA.  Buot,  “  Operation  Principle  of  Resonant  Tunneling  THz 
Oscillator  at  Fixed  Bias  Voltages”,  to  appear  in  VLSI  Design,  2001. 

3.  P.  Zhao,  H.L.  Cui,  D.L.  Woolard,  F.A.  Buot,  and  K.  Jensen,  “  Self-induced,  sustained  electron  oscillation 
of  a  resonant  tunneling  diode”,  submitted  to  Physical  Review  Letters,  (2001). 

4.  H.  J.  Quan,  X.M.  Weng,  and  H.L.  Cui,  “Effects  of  Anisotropy  on  the  Electron  Transport  in  4H-SiC”, 
Physica  Status  Soilidi  (b)  219, 339-346  (2000). 

5.  P.  Zhao,  H.L.  Cui,  D.L.  Woolard,  K.L.  Jensen,  and  F.A.  Buot,  “Origin  of  hysteresis  and  plateau-like 
behavior  of  the  I-V  characteristics  of  resonant  tunneling  diodes”.  International  Journal  of  Modem  Physics 
B,  14,411-426(2000). 

6.  P.  Zhao,  H.L.  Cui,  D.  Woolard,  K.  L.  Jensen,  and  F.  A.  Buot,  “Simulation  of  resonant  tunneling  structures: 
Origin  of  the  I-V  hysteresis  and  plateau-like  structure”,  Journal  of  Applied  Physics,  87,  1337-1349(2000). 

7.  F.A.  Buot,  P.  Zhao,  H.L.  Cui,  D.  L.  Woolard,  K.  L.  Jensen  and  C.M.  Krowne,  “Emitter  quantization  and 
double  hysteresis  in  resonant-tunneling  structures:  A  nonlinear  model  of  charge  oscillation  and  current 
bistability”,  Physical  Review  B  61,  5644-5665  (2000). 

8.  P.  Zhao,  H.L.  Cui  and  N.J.  Horing,  “Quantum  Transport  Equations  for  single-band  and  Multi-band 
Semiconductors”, Philosophical  Magazine,  B80. 1359-1392,  (2000). 

9.  P  Zhao,  HI.  Cui  and  D.L.  Woolard,  “Numerical  simulation  of  G-V,  C-V ,  and  L-V  characteristics  of  a 
resonant  tunneling  diode”.  International  Semiconductor  Device  Research  Symposium,  December,  1999, 
Charlottesville,  VA. 

10.  P.  Zhao  and  H.L.  Cui,  “Quantum  transport  equation  for  a  two-band  system”,  Physics  Letters,  252, 243-247 
(1999). 


11.  X.M.  Weng,  H.L.  Cui  and  HJ.  Quan  “Simplification  of  scattering  rates  calculation  in  nonparabolic 
hydrodynamic  equations”.  Journal  of  Applied  Physics,  85, 2658-2662  (1999). 

12.  C.C.  Lee,  H.L.  Cui,  and  J.  Cai,  “Transient  simulation  of  a  Si-MESFET  using  hydrodynamic  balance 
equations”.  International  Symposium  on  VLSI  Technology,  Systems,  and  Applications,  1998. 

13.  C.C.  Lee,  H.L.  Cui,  J.  Cai,  R.  Pastore,  D.  Woolard,  and  D.  Rhodes,  “Two-dimensional  hydrodynamic 
device  modeling”.  Symposium  on  Multi-Simulations,  1998. 

14.  X.L.  Lei  and  H.L.  Cui,  “Balance  equations  for  electron  transport  in  an  arbitrary  energy  band  driven  by  an 
intense  terahertz  field  -  application  to  superlattice  miniband  transport”,  European  Physics  Journal  B,  4, 513 
(1998). 

15.  J.  Cai,  H.L.  Cui,  E.  Lenzing,  R.  Pastore,  D.  Rhodes,  and  B.  Perlman  “Hydrodynamic  device  modeling  with 
band  nonparabolicity”,  VLSI  Design,  6, 181-184  (1998). 

16.  X.L.  Lei,  N.J.  Horing,  and  H.L.  Cui,  “Microscopic  three-dimensional  theory  of  space-charge-wave  modes 
in  GaAs-based  planar  superlattices”,  Communications  of  Theoretical  Physics,  29,  515-522  (1998). 

17.  X.M.  Weng  and  H.L.  Cui,  “Simulation  of  a  submicron  3C-SiC  ballistic  diode  based  on  the  Lei-Ting 
hydrodynamic  balance  equations”.  Journal  of  Applied  Physics,  83, 3129-3133  (1998). 

1 8.  X.M.  Weng  and  H.L.  Cui,  “Numerical  simulation  of  3C-SiC  ballistic  diode  based  on  the  Lei-Ting 
hydrodynamic  balance  equations”.  Proceedings  of  International  Semiconductor  Device  Research 
Symposium,  December,  1997,  Charlottesville,  Virginia. 

19.  X.L.  Lei  and  H.L.  Cui,  “Suppression  of  absolute  instability  by  elastic  scattering  in  semiconductor 
superlattices”,  Journal  of  Physics:  Condensed  Matters,  9, 4853-4861  (1997). 

20.  X.M.  Weng  and  H.L.  Cui,  “Temperature  dependence  of  Hall  mobility  of  SiC”,  Journal  of  Physics: 
Condensed  Matter,  9, 7089-7094  (1997). 

21.  X.M.  Weng  and  H.L.  Cui,  “High  temperature  hot  electron  transport  in  6H-  and  3C-SiC”,  Physica  Status 

~  Solidi  (b),  201, 161-166  (1997). 

22.  X.L.  Lei,  H.L.  Cui,  N.  Horing,  and  K.K.  Thomber ,  “Superlattice  space-charge-wave  analyses  in  a 
modified  drift-diffusion  model  and  in  a  simplified  balance-equation  formulation”,  Z.  Phys.  B  104, 221-225 
(1997). 


